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^ ■ Abstract 



A general non-perturbative analysis of the renormalization properties of A/ = 3/2 
four-ferniion operators in the framework of lattice regularization with Wilson fermions 
a I is presented. We discuss the non-perturbative determination of the operator renormal- 

ization constants in the lattice Regularization Independent (RI or MOM) scheme. We 
also discuss the determination of the finite lattice subtraction coefficients from Ward 
Identities. We prove that, at large external virtualities, the determination of the lattice 
mixing coefficients, obtained using the RI renormalization scheme, is equivalent to that 
based on Ward Identities, in the continuum and chiral limits. As a feasibility study of 
our method, we compute the mixing matrix at several renormalization scales, for three 
values of the lattice coupling /?, using the Wilson and tree-level improved SW-Clover 
actions. 



1 Introduction 

The renormalization constants of lattice operators are necessary ingredients in the 
prediction of physical amplitudes from lattice matrix elements. Schematically, the 
physical amplitude A^^p associated to the transition a ^ P induced by the composite 
operator O is given, in the formalism of the Operator Product Expansion (OPE), by 

A^^p = Cw{^^)Zo{a^^){a\0{a)\P) (1) 

where Cw is the Wilson coefficient, /x is the renormalization scale, a is the lattice 
spacing, {a\0{a)\P) is the matrix element of the bare lattice operator relevant to the 
physical process and Zq its renormalization constant]^. The Wilson coefficient Cvi/(/i) 
is calculated in perturbation theory at the renormalization scale /x. It depends on 
the renormalization scheme (a mass-independent scheme is implied throughout this 
paper). Wilson coefficients are known to NLO; see refs. |]^-|^]. The matrix element 
must be calculated non-perturbatively; the only known method for computing it from 
first principles (at a fixed cutoff a~^) is lattice QCD Q. The operator renormalization 
constant Zq is the link between the matrix element, regularized on the lattice, and its 
renormalized, continuum counterpart. So far, three methods have been implemented 
to calculate it: 

1. Lattice Perturbation Theory (PT) §- g. 

2. Lattice Ward Identities (WI), in the spirit of ref. ]T0| , 11], either with hadronic 
states (see |12|-[0]), or with quark states (see |]l4|,]l5| and, more recently []l8|). 



or even with the Schrodinger functional (see |1^). The WI method is non- 
perturbative. 

3. The Non-Perturbative (NP) method of |2C], which consists in imposing renor- 



malization conditions at quark correlation functions with operator insertions, at 
a given scale (see Q-Q). 



Once the systematic errors of these calculations are understood and kept under con- 
trol, the matrix elements, computed in simulations, can be properly renormalized and 
reliable non-perturbative estimates of physical quantities can be obtained. 

The simplest case is that of the matrix elements of two-fermion (dimension-three) 
operators, from which, for example, the decay constants of light vector and pseudoscalar 
mesons can be extracted. Several complications are avoided due to the fact that the 
operators renormalize multiplicatively. The situation is less straightforward in the 
case of four-fermion (dimension-six) operators. Their renormalization is also additive; 
mixing occurs with other operators of the same dimension and, in some cases, also 
with operators of lower dimensionality. These operators must then be subtracted in 
order to make the original operator finite. Mixing with lower-dimensional operators 



^ Here our notation is schematic: in general we have a set of operators and Zq is a matrix which mixes 
those operators which form a complete basis under renormalization. 



(e.g in the problem of understanding the AI = 1/2 rule) is characterized by several 
theoretical and technical subtleties, which have been addressed in [25|. 



This paper consists in a study of the renormalization properties of the complete 
basis of dimension-six four-fermion operators in the absence of power subtractions 
(i.e. operators mixing only with others of equal dimension), using non-perturbative 
methods. We will discuss below several phenomenological applications of this problem, 
related to light (strange) flavour phenomenology. We point out however, that since we 
are dealing with a mass-independent renormalization scheme, our results can also be 
applied to heavy-flavour phenomenology (e.g. the A.B = 2 process). This calculation 
is also of relevance to the evaluation of the A/ = 1/2 matrix elements, since the 
anomalous dimension of the octet operators is unchanged by the mixing with lower 
dimensional operators (and the related power subtractions) [p^ ]. The most common 
examples, concerning ii'-meson decays, which involve operator mixing without power 
subtractions are the following (note that SU(2)-isospin breaking effects are neglected): 

• The study of K^-K^ oscillations involves the computation of Bk (the B-parameter 
of the i('-meson), which is obtained from the AS* = 2 matrix element {K^\0 '^\K^) 
of the operator 

o^^=2 = {sjf:d)is^f:d) (2) 

where s and d stand for strange and down quarks. Our conventions for Dirac and 
colour matrices, indices etc. are given in Appendix ^. 

• The study of the A/ = 3/2 sector of the decay K — > vrvr involves the computation 

3/2 

of the matrix elements (vrvrlOg ^^QJ-fC) of the operators 

or = {s^ifid^){u''iX) + {-s^iX){u^jj:d^) - {s^7fid^){d^jf:d^) 

0%' = i-s^ll:d^){n^jl:n^) + (.-^7>^)(n^7,^d^) - (^-^7,^d^)(d"^7,^d^) (3) 

with u standing for the up quark and A, B denoting colour indices. The above 
operators transform as members of the (27, 1) representation of the SU{3)l ® 
SU{3)ji chiral group. Moreover, the study of the electropenguin contribution to 
the decay K —^ vrvr involves the computation of the matrix element (vrvrlOy g \K) 
of the operators 

O?/^ = {s^lfld^){u^jX) + («-^>^)(^^.^^^) - i^^lfld^)id''l^d^) 

o'f = {-s^il:d^){n^^X) + («-^7>^)(^^7,^d^) - i-s^il:d^)id^7X) (4) 

which are also AI = 3/2, but transform as an (8, 8) representation of the chiral 
group. These K — > vrvr matrix elements are essential to the calculation of e'/e. 
They can be obtained, through soft pion theorems, from the matrix elements 
(vr^^lO^ 1^^) (for k = 7,8,9,10). From the single-state matrix elements, the 

3/2 

corresponding i?-parameters, i?y g , can be extracted. Note that, in the limit of 



3/2 

degenerate quark masses, both operators OgiQ, having a "left-left" Dirac struc- 
ture, renormalize in the same way as the O ^^ operator. Moreover, the matrix 
elements (vr+IOg^Ql-ftT"'") have the same i?-parameter as {K^\0 "^{K^) . These 
matrix elements should vanish in the chiral limit. On the contrary, the operators 

3/2 

O7 g , having a "left -right" Dirac structure, obey different renormalization prop- 
erties. The corresponding matrix elements (vr+IO^g \K^) do not vanish in the 
chiral limit 

Important information on the physics beyond the Standard Model, such as SUSY, 
can be obtained by studying FCNC processes and, in particular, AF = 2 transi- 



tions (see [27 1 and references therein for a discussion). Besides the i?-parameters 
of the operators listed above, such processes also require the knowledge of the 
i?-parameters of the operators 

(s-^(l-75)d^)(s^(l-75)d^) 

(s-^(l-75)d^)(s-^(l-75)d^) 

(s-^(l-75)rf^)(s-^(l + 75)d^) (5) 

(s-^(l-75)d^)(s^(l+75)d^) 

Recent lattice results on all these S-parameters (with Wilson fermions) can be found 
inrefs. ||, ||, ||, |0). 

The matrix elements discussed above are extracted from the large time asymp- 
totic behaviour of three-point correlation functions of the form {P^O '^Pk) and 
{PJ^Oj Ptt) (for k = 7,8, 9, 10) where Pk and P^^ are pseudoscalar sources of suitable 
quark flavour. Expressed in terms of traces of quark propagators, these correlation 
functions correspond to the so-called "eight"- shaped quark diagrams. Any "eye"- 
shaped contributions of these correlation functions cancel in the limit of degenerate up 
and down quarks, which is the case under study. 

The paper is organized as follows: in sec. ^ we illustrate, as a concrete example of 
the problems arising in lattice renormalization, the mixing of the operator O . Its 

matrix element {K^\0 "^IK^) is very sensitive to the various systematic errors which 
affect its chiral behaviour and therefore provides a good case-study. In sec. [^ we discuss 
in full detail the problem of mixing under renormalization of all four-fermion dimension- 
six operators. The operators are classified according to their discrete symmetries and 
an operator basis, convenient for our purposes, is chosen. Subsequently, the operator 
mixing under renormalization is divided into two parts: mixing that would occur even 
if all the symmetries were respected by the regulator, and mixing which is induced 
on the lattice by the chiral symmetry breaking of the Wilson term. In sec. ^ the NP 
renormalization is applied to the cases at hand. The renormalization conditions of the 
so-called RI (or MOM) scheme are worked out in terms of projected amputated Green 
functions, for the complete operator basisQ In sec. |5| we derive WI's suitable for the 



^For continuum calculations in the RI scheme and the anomalous dimensions of the complete basis of 
operators, calculated in PT to NLO, see ref. g]. 



determination of the lattice mixing coefficients in the chiral hmit. We show that, in 
this respect, the RI and the WI methods are equivalent in the region of large external 
momenta. We also show the independence of the lattice mixing coefficients from the 
renormalization scale. Finally, in sec. ^ we present the renormalization constants and 
lattice mixing coefficients, obtained with the NP method, with both the Wilson and the 
tree-level improved Clover action at several values of the lattice bare coupling /3. Some 
of these results have been used in refs. [^, ^] in order to derive fully non-perturbative 
estimates of the 5-parameters. 

2 Renormalization of 0^^='^ and systematic er- 
rors 

The testing ground for the restoration of chiral symmetry in the continuum limit of lat- 
tice computations has been the chiral behaviour of the matrix element {K^\0 '^\K^) , 
which, if properly renormalized, vanishes when the i^-meson becomes massless |31|| . 
Although early attempts with Wilson fermions |32]-[^] have given reasonable mea- 
surements of Bk, they have been less successful in controlling this chiral behaviouiQ. 
The root of the problem is the operator subtraction outlined above, which we now 
discuss in some detail: O "^ mixes with other operators Oj of the same dimension 
but with "wrong naive chirality". Thus, the (;U-dependent) K^-K^ matrix element of 
the renormalized operator O is given in terms of the (a-dependent) bare matrix 

elements by: 






0^'=\a)+y^Zt'=\9l)0,{a) 



\K^) 



(6) 

The overall renormalization constant Zq "^{a^, Qq) is logarithmically divergent whereas 
the Zj^'^^^(5rg)'s are finite mixing coefficients; (7o(a) is the lattice bare coupling (also 
expressed as /? = Q/qq)- All Z's can be calculated, at least in principle, in perturbation 
theory (PT), provided that /i, a~^ ^ Aqcd- The one-loop perturbative calculation of 
the Z's, both for the Wilson and Clover lattice actions, has been performed in [^- Q. 
Although the renormalized operator O (/i) is constructed so as to have the correct 

chiral properties in the continuum limit, the implementation of operator mixing in 



early works [^2[-[34] was subject to two main systematic errors: 

1. The overall renormalization constant Zq^'^ and the mixing coefficients Z^^"^ 
were calculated in one-loop PT. Thus, they suffered from ©(^g) systematic errors. 



^Bk has also been obtained with staggered fermions |35| mainly in the quenched approximation (see 
ref. Ipql for a review). The (surviving) chiral symmetry in the staggered fermion formalism ensures the 
vanishing of the relevant matrix element in the chiral limit. 



2. The bare matrix elements (K°|0^'^-2(a)|K°) and {K^\Oi{a)\K°) were calculated 
non-perturbatively (numerically) with the Wilson action at fixed cutoff a (i.e. 
fixed coupling (^q). They were therefore subject to 0{a) systematic errors. 

Both sources of systematic error may be considered responsible for the non-vanishing 
of the matrix element {K^\0 '^\K^) in the chiral limit. In order to reduce these 
errors, the following remedies have been proposed: 

1. Boosted Perturbation Theory (BPT) should, according to [^], improve the be- 
haviour of the perturbative series of the renormalization constants. The system- 
atic error due to the 1-loop truncation of the BPT estimates of the Z's remains 
C(5o)> t>u^ it is hoped that it is smaller than the C'(g'o) error of the original PT 
expansion. 

2. Discretization errors can be reduced by using improved actions and operators. 
To this purpose, besides the Wilson action, the Clover action |^^ has also been 
implemented in the calculation of the weak matrix elements {K^\0 "^ {a)\K^) 
and {K^\Oi{a)\K^). Matrix elements of tree-level improved operators, calculated 
with the tree-level improved Clover action, have discretization errors of 0{gQa), 
instead of 0{a) as in the Wilson casd^. 



3. Using the non-perturbative (NP) method of [^, 21|, a more accurate evaluation of 
the renormalization constants can be achieved. The method consists in imposing 
renormalization conditions directly on quark four-point Green functions with op- 
erator insertions. Higher-order contributions, including non-perturbative effects, 
are taken into account. With this method, the Z's have the same discretization 
errors (due to the finiteness of the lattice spacing) as the operator matrix elements 
we are ultimately interested in. Another non-perturbative method consists in the 
determination of the mixing coefficients from WI's on quark states ||l^. We stress 
that the WI method can only be applied to the lattice mixing coefficients; the 
overall renormalization Zg "^^^(a//, ^q) can only be determined non-perturbatively 
from NP methods; e.g. those of refs. |l^, |2^, |2^. Both WI and NP results are 
expected to have similar discretization errors. 

The systematic errors of the operator renormalization described by eq. (^) depend upon 
the method of calculation of its various terms: if PT (or BPT) is used for the Z's and 
the Wilson action for the matrix elements, the renormalization is good to O^Qq) and 
0{a). If PT (or BPT) is used for the Z's and the tree-level Clover action is used for 
the matrix elements, the renormalization is good to ©(^q) and 0{gQa). If NP (or WI) 
renormalization is used for the Z's then the error is 0{a) or 0{gQa), depending on the 
action used (Wilson or tree-level Clover). Of course, the same considerations are valid 
for the renormalization of the operators of eqs. (p|)"(|5|). 



* Recently, it was shown that the non-perturbative Clover improvement, proposed in p9| reduces the 
above error to 0{a^). We have not implemented this version of the Clover action in the present work, 
because it involves more complicated operator mixing. 



A study of the effect of each of these three improvements had first been made 
in ref. [^], in which a direct comparison of the various errors in the computation 
of the AS = 2 matrix element was carried out. The main resuh of [^9[ is that the 
chiral behaviour of the renormahzed {K^\0 '^\K^) is not improved significantly if 
BPT is used for the calculation of the mixing constants Z^^"^ and the Clover action 
is implemented in the computation of the matrix elements^. What really makes a 
difference is the NP computation of the Z's, which was shown to improve significantly 
the chiral behaviour of the matrix element. This conclusion has been confirmed in 



18 1 with Wilson fermions and for several f3 values. We stress that in |l^] the non- 



perturbative computation of the mixing Z's was performed using WI's on quark states, 



as opposed to the NP method of [|39|. 



In all early works, ref. |3^ included, the renormalization of O consisted in 

the mixing with only three dimension-six operators with the "wrong naive chirality" . 
This mixing had been explicitly derived in Q, using lattice PT at one- loop. However, 
by applying the general symmetry arguments of |41| to the case in hand, one can 



show that a fourth operator should be added in the sum of eq. (y), as pointed out 
in refs. |18, 24, 33]. This operator, absent at the one-loop level of PT, has a mixing 
constant which is at least O((7o)- ^^ gives a small (but non-negligible) contribution 
to the restoration of the chiral behaviour of the matrix element ||2^. In refs. |24, 29 1 
and the present work, this operator, omitted in ref. |39], has been properly taken into 
account. 



3 Mixing of dimension-six four-fermion opera- 
tors on the lattice 

In this section, we introduce the basic notation and study the mixing of generic four- 
fermion operators in the presence of explicit chiral symmetry breaking, induced by the 
Wilson term. This implies that mixing with operators with the "wrong naive chirality" 
is allowed, even in the chiral limit. Following |^1| (see also ||3^ ), we obtain the complete 
basis of dimension-six, four-fermion operators which mix under renormalization, relying 
on general symmetry arguments based on the vector flavour symmetry, which survives 
on the lattice. In order to solve this problem, it is convenient to work with massless 
fermions with four distinct flavours tpj^ / = 1, . . . , 4. The correlation functions arising 
from these operators do not contain "eye" -shaped contributions and therefore their 
renormalization properties are identical, up to trivial factors, to those of the physical 
AS" = 2 and A/ = 3/2 operators considered in sect. |l| (recall that we are assuming 
SU(2)-isospin symmetry). The introduction of "light" masses (i.e. masses which are 



much smaller than the renormalization scale) will be discussed in subsect. 5^. We will 
not address the problem of "heavy" masses (of the order of, or greater than the typical 
scales and the UV cutoff) in this work. 



For a different conclusion see ref. WOl 
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Table 1: Classification of four- fermion operators Or{i)r{2) according to V and useful products 
of their discrete symmetries C, S' and S" . These properties are also valid for the operators 
Ojar(i)tar{2). For the operators Op{i)p(2) 
columns CS' ^^ CS" and CVS' ^^ CVS" 



and 0^p(i)^^p(2), we must exchange the entries of the 



Note that Off 



Ott and O 



TT 



o. 



TT- 



We define the generic four fermion operators 






(V5it''r(^)^4)(^3i"r(2)v;2 



(7) 



where r'^^ and T'^-* denote any Dirac matrices and i" are the colour-group genera- 
tors (referred to as colour matrices in this work). For notation and conventions, see 
Appendix ^. 



3.1 Operator classification according to discrete symme- 
tries 

Let us start by addressing the problem of mixing of the generic operators of eq. (|^ 
under renormalization. First of all we stress that they cannot mix with operators of 
lower dimensionality, because such operators do not have the four-flavour content of the 



original ones^. Thus, Op(i)p(2) can mix with any other dimension-six operator, provided 
it has the same quantum numbers (i.e. with any operator which has the symmetries of 
Or(i)r(2) and of the action). The generic QCD Wilson lattice action with 4 degenerate 
quarks is symmetric under parity V, and charge conjugation C. Moreover, there are 
three other useful (flavour) symmetries of the action, namely the flavour exchange 
symmetry iS = (^"2 ^^ V'4) and the switching symmetries S' = {ipi <-^ ^2) "03 ^^ ^^4) 
and S" = [ipi ^^ 04, V'2 <-> "03) @l- III Table || we classify the operators Or(i)r{2), for 
all r^^) and T^^) combinations of interest, according to the discrete symmetries V, C, 
S' and S" . Parity violating operators (middle section of the Table), for which CS" is 
not a symmetry, have been (anti)symmetrized in order to obtain eigenstates of CS" 
(bottom section of the Table). We adopt the notation 

0[r{i)r(2)±r(2)r(i)] = O-pwrd) ± Op(2)r{i) (8) 

Note that the results of Table || apply also to the operators 0^ap{i)^ap{2) since, upon 
performing the symmetry transformations, sign differences, resulting from the presence 
of the colour t"" matrix, disappear because the colour matrices appear quadratically. On 
the other hand, O^^j (2) is obtained by applying S on Op(i)p{2) . Since S transforms S' 
into S", the properties of Table |l| also apply to O^^j (2) ; but with all S' and S" columns 
exchanged. Again, the operator 0^p(i),ap(2) has the same properties as O^^j (2)) since 
the colour matrix t" appears quadratically. 

Our aim is to find complete bases of operators which mix under renormalization. 
Thus, the first task is the elimination of the operators which are not independent. In 
particular, given the mixing of Op{i)p(2) with all other allowed Op(i)r(i) 's and Op{i)p(j) 's, 
we do not need to consider also the mixing with the OtT(»)t°r(j)'s and 0£^„(j) „(j)'s, 
since they can be expressed in terms of the Op(i)p(j)'s and the O^ijpj^j's. This can 
be seen by applying the standard identity of colour matrices (Fierz transformation of 
colour indices) 

tABtcD = -l^^AB^CD + -j^^AD^CB (9) 

on the t'^'s of a given operator. For the operator O^apCij^apO) the result has the general 
form 

Ojap(i)tap(j) = — ■^^Or(»)r(j) + :^ 2^ C^„mOp(„)p(„) (10) 

c n,m 

where the sum runs over all the Dirac matrices obtained by the Fierz transformation 
of r'^^r^^ The factors Cnm are the appropriate constants of the Fierz transformation 
of the Dirac matrices (see Appendix y for details on Fierz transformations in Dirac 
space). Analogously we can express 0^p{i)jap{j) in terms of Op(i)pO) and Op(i)p{j)- 
Therefore, in the following, it is adequate to limit ourselves to the mixing of Op(0r(j)'s 
and Op(j) yj's, which form a complete basis of 20 independent operators (these are the 



^This statement remains true for the operators of eqs.(||)-(|^). The A/ — 1/2 operator however, mixes in 
general with operators of dimension- five and -three. 



5 operators of the top section of Table Q, the 5 operators of the bottom section and 
their 10 Fierz counterparts). 

Furthermore, by classifying the operators according to the discrete symmetries 
listed in Table 0, the original basis of 20 operators can be further decomposed into 
smaller independent bases. An immediate decomposition is that into two bases, of 10 
operators each, with definite parity ("P = ±1). Further decompositions occur upon 
using the remaining CVS' , CVS" and S symmetries. These we now perform case by 
case. 

3.2 Parity conserving operators 

All of the parity conserving operators Orr are eigenstates of all the discrete symmetries 
listed above, with eigenvalue +1. Therefore, each of them can mix with the other four, 
and also with the five Opp's; the complete basis consists of 10 operators. We now 
rotate our basis into a new one, consisting of 10 operators which are also eigenstates 
of S with eigenvalues ±1: 



Otr 



1 



1 



Opp ± ofpj = - [(V^ir^2)(V^3rV4) ± (V5irv4)(V53rV2)] (n) 



Clearly, the 5 Opp's, corresponding to 5 = +1, mix only among themselves; the same 
is true for the Opp's which have S = —1. Thus the original basis of 10 operators has 
been decomposed into two independent bases of 5 operators each. 

In this work, since we are eventually interested in the renormalization of the oper- 
ators of eqs. (|2|)~® ) we opt for the basis 



Qf 


- ^[VV+AA] 


Qt 


- ^[VV-AA] 


Qt 


- '^[SS-PP] 


Qi 


- '^[SS+PP] 


Qf 


— Urpj, 



(12) 



Other choices of basis can be found in [g, |l^] . All such bases are simply linear combi- 
nations of each other. 

3.3 Parity violating operators 

The parity violating four- fermion operators listed at the bottom section of Table |l| 
do not all have identical CVS' and CVS" values. We will now establish their mixing 
pattern, with the aid of S symmetry. 

The operator 0\yA+AV] mixes with 0^^_(_^yi only, forming a basis of two operators 
with CVS' = CVS" = +1. In analogy to the parity conserving case, we rotate this 



basis into 



'^[VA+AV] 



o 



[VA+AV] 



±0 



[VA+AV] 



(13) 



Since O^^ , ^y, has 5 



+1 and Or^^ , ^yi has 5 = —1, they do not mix with each 
other; the two operators renormahze multiphcatively. In other words, the original basis 
of two operators has been decomposed into two bases of one operator each. 

The operators 0[yA-Ay] ^-^id ^\sp-ps] ™^^' ^ii^ce they both have CVS' = +1 and 
CVS" = — 1. Similarly, Oiy^.^yi and 0[sp_ps] form another basis, with CVS' - 
CVS" = +1. It is convenient to combine these two bases into 4 operators: 



1, 



^[VA-AV] 
'^[SP-PS] 



1 r 

2 
1 

2 



(0[ 



VA-AV] 



±o 



[VA-AV] 



0[sp-ps] ± Oygp_pgT^ 



(14) 



None of these operators have definite CVS' or CVS" . However, they have definite 
S = ±1. There is mixing between the two 5 = +1 operators (0|y^_yiv'l' ^tsP-PS") 
and also between the two S = —\ operators (0^^_^yi, ^rsp-psi)- 

We follow a similar line of reasoning for the operators 0[5P4.p5], Ofgp_^pg^, Orpf 
and O^^ (all have CVS' = CVS" = — 1). We rotate the basis of these four operators 



TT 



into 



^[SP+PS] 



-1± 



Oh 



0[SP+PS] ± 0[_5.p^p5.] 



Oj.f ± O. 



TT 



(15) 



Thus, once more, the original basis of the four operators of eq.(O) can be decomposed 
into two bases, of two operators each, with definite 5 = ±1. 

For reasons of notational compactness, we will use the following redefinitions: 






Q-. 



Qt 



Q 



^[VA+AV] 
^[VA-AV] 

^[SP-PS] 
^[SP+PS] 



(16) 



3.4 Operator subtraction and chiral symmetry 

This completes the discussion of the mixing of dimension-six operators on the lattice, 
in the general case of 4 distinct flavours (degenerate or not). We now summarize the 
result to keep in mind. For the parity conserving operators, the renormalization is 
given by 

Qf = Zt^Qf (i,j = l,...,5) (17) 
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whereas for the parity violating ones we have 



Q 






± 



Vi',3 



,5) 



(18) 



where Qi aiid Qi 



are the bases of renormaUzed operators and Zi 



renormaUzation matrices (summation over repeated indices is imphed). 



and Z;^ are the 



± • 



On the basis of CVS symmetries, it is important to notice that the matrix Z^- is 



«j 



a (relatively sparse) block diagonal matrix. We show this explicitly by re-writing the 
last equation according to the allowed mixing of the operators derived in the previous 
Subsection: 



( Qi\ 

Q3 

Q4 

VQ5/ 



( Zii \ 


HI 


( 2i \ 


Z22 ^23 




Q2 


^32 ^33 




Q3 


2:44 2:45 




Q4 


iy Z54 ^55 / 




^25^ 



(19) 



The abbreviated notation with the ± superscript of the column vectors and matrix 
should be transparent to the reader. In conclusion, the lattice does not induce ex- 
tra subtractions for the parity violating sector (Q^ ; k = 1, ... ,5), since the mixing 



occurring according to the pattern of eq. (|19D, is also valid in the x^S scheme; see 
below. 

In the case of parity conserving operators, we proceed in a different way. We find 
it convenient to separate the above operator mixing into two classes: 1) the first is the 
lattice subtraction, which consists in correcting the operator mixing, induced by the 
breaking of chiral symmetry due to the Wilson term of the action; 2) the second is the 
subtraction which survives in the continuum limit. In order to facilitate this separation, 
let us suppose that there is a regularization scheme which, unlike the lattice, respects 
chiral symmetry; i.e. the regularized action has no chiral symmetry breaking term 
(recall we are working at zero quark mass). We will call this hypothetical scheme the 
Chirally-symmetric Regularization Scheme {xRS for short). In this scheme, we can use 
chiral symmetry in order to establish some extra selection rules for the renormalization 
of the operators of interest. Although the resulting "operator renormalization" is 
incomplete on the lattice, it is the one which should be recovered in the the continuum 
limit, since the renormalized theory must have the desired chiral properties. The 
remaining lattice subtractions are due to the presence of the Wilson term. 

In order to derive the selection rules in the x^S scheme, it is adequate to consider 
two discrete axial symmetries. The first (denoted by X24) acts only on flavours 2 and 
4: 



V'2 
V'4 



Z75V'2 



1p2 



ilp2l5 



(20) 
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Symmetry 


Qt 


Qt 


Qt 


Qt 


Qt 


Qr 


Q2 


Q3~ 


Ql 


Qi 


X24 


-1 


+1 


+1 


-1 


-1 


-1 


+1 


+1 


-1 


-1 


Xl2 


+1 


Q2 


-Q3 


-1 


-1 


+1 


Q^ 


-Qt 


-1 


-1 



Table 2: Classification of four-fermion parity conserving operators Qf^ (A; = 1, . . . , 5), accord- 



ing to the discrete symmetries X24 and Xi2- The parity violating operators Q^ {k 
transform in the same way. 



.5) 



The second symmetry (denoted as X12) acts only on flavours 1 and 2: 



■01 
V'2 



Z75V'2 



1p2 



#175 



(21) 



Under these transformations, the four-fermion operators transform as shown in Table |2|. 
The symmetry X12 implies that Q^ renormalize multiplicatively, whereas QJ and Qt 
{Q^ and Q^) mix with each other. From the X24 symmetry we deduce that Qt and 
Qt iQ2 ^^*^ Q3) ^1^0 ^^ with each other. In the absence of explicit chiral symmetry 
breaking, we conclude that the mixing structure is the same as that considered above 
for the parity-violating counterparts. The corresponding parity conserving operators 
{Qj^ ; A: = 1, ... ,5), belonging to the same chiral representations, would have a xRS 
mixing pattern similar to that of eq. (|l9|): 



f Qi\ 
Q2 
Qs 

Qa 

\Q5 J 



f Zii \ 

Z22 Z23 

Z32 Z33 

Z44 Z45 

V Z54 Z55 j 



( Qi\ 

Q2 

Q3 

Q4 
V Q5/ 



(22) 



where the Qi represent the bare operators in the x^S scheme. In the presence of the 
Wilson term, the Qi are the lattice subtracted operators defined as 



Q2 

Qs 

Qa 

\Q5 J 



f Qi\ 
Q2 
Q3 
Qa 

\Q5 J 



+ 



( ° 


A12 


Ai3 


A 14 


A21 








A24 


A31 








A34 


A41 


A42 


A43 





I A51 


A52 


A53 






Al5 \ 

A25 
A35 


0/ 



± 



f Qi\ 
Q2 
Q3 

Qa 



(23) 



In other words, first the lattice subtraction is performed, followed by the x^S renor- 
malization. The above mixing pattern is abbreviated, in matrix form, as 



Q^ 
Q^ 



Z^Q^ 

[/ + A±]Q 



± 



(24) 
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where / is the 5x5 unit matrix and the subscript x stands for x^S subtractions \\. 
The renormahzation of the parity conserving sector is given by Z = Z^[I + A ]. 
Note that, using continuous chiral transformations, in the hypothetical x^S, it is easy 
to show that A^ = and Zij = Zij. 

4 NP renormalization: the lattice RI renormal- 
ization scheme 

In this section, we outhne the strategy used in the determination of the renormahzation 
constants, and discuss the appropriate renormahzation conditions. According to the 
NP method proposed in [^, the renormahzation conditions are imposed in momentum 
space on the projected amputated Green's functions. In this work we will always 
consider Green's functions where all external quark legs have the same momentum p. 
This choice is, of course, not unique but is the simplest way to regulate the infrared 
divergences. An important point is that the renormalization conditions need to be 
imposed at large Euclidean p"^. In the determination of the mixing coefficients this 
condition enables us to neglect soft or spontaneous chiral symmetry breaking effects 
which are not induced by the Wilson term and cannot be computed. As demonstrated 
in sec. |5|, only at large external momenta is the NP method, described below, equivalent 
to the WI one. The large p'^ condition is also necessary in the evaluation of the 
overall renormalization constants, since the standard procedure of obtaining physical 
amplitudes from renormalized matrix elements requires perturbative matching with 
the Wilson coefficients in the continuum, at fixed gauge coupling. We will describe in 



subsect. 4.3 an alternative procedure, for which this condition can be relaxed. 

The renormalization method presented in this section is an extension of the one 
introduced in [21| for the renormalization of the O four-fermion operator. As in 



ref. [£^], we will define projected-amputated Green functions, on which suitable renor- 
malization conditions will be imposed. With respect to ref. [21 1, the novelties consist 



in using a complete operator basis for the parity-conserving sector and in extending 
the method to the parity-violating operators. 

4.1 Amputated Green functions and their projectors 

We first give some general definitions. Since the non-perturbative renormalization 
conditions are to be imposed on quark states, we shall need the general expression 
of the four-point Green's functions of the operators O (^j (j)- Denoting by xi, 2:3 and 
3:2, 2:4 the coordinates of the outgoing and incoming quarks respectively, we define the 



^Attention is drawn to the notation adopted from now on in this paper: whereas the elements of matrices 
Z^ and A^ are denoted, in standard fashion, as Zfj and A^ , those of matrix Z^ are denoted by Z^ (i.e. 
the subscript x is dropped for notational economy). The elements of matrix Z^, appearing in eq. (|l7|), will 
never be used from now on. 
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connected one-particle-irreducible Green's functions as 

G^(1)P(2){X1,X2,X3,X4) = (V'l(xi)^2(2;2)Op(i)p(2)(0)V'3(2;3)V'4(a;4)) , (25) 

where (• • •) denotes the vacuum expectation value, i.e. the average over the gauge-field 
configurations. The generic four-fermion operator Op(i)p{2) , placed at the origin, is 
written as 

O^(i)r(2)(0)= (26) 

Given the complexity of the notation, we shall elaborate on this formula. In the case 
of parity conserving operators, the two Dirac matrices are equal; i.e. we are dealing 
with the Opp of eq. ([TT|). The case of two different Dirac matrices F'^' ^ r*-^-* applies 
to the parity violating operators of eqs. (p^) , (p^) and (jl^), for which, for example, 
r(^) = V and r*-^) = A. Note that G (^j (2) depends implicitly on the four colour and 
Dirac indices carried by the external fermion fields. These will be shown explicitly 
from now on. Following the conventions of Appendix |A|, upper-case Latin superscripts 
{A, B, . . . , R, S, . . .) will denote colour indices of the fundamental SU(3) representation. 
Greek lower-case subscripts (a, /?,..., p, cr, .. .) will denote spinor indices. The letters 
m and n are reserved for the Lorentz indices of the Dirac matrices. Repeated Lorentz 
indices are summed according to the convention adopted in Appendix ^. 

The Fourier transform of the non-amputated Green's function (|25|), at equal exter- 
nal momenta p, has the form 

G^,.,,,.Ma^,T = \ [(r«-(p)i^F(2)-(p),T) T (F(^)-(p);^f F(^)-(p),^/)] , (27) 

where 

r«"^(p)^| = S{mif^fr[l^S{pWl^\^l {^ = 1, 2) (28) 

and S{p\0) is defined by 

S{p\0)= f d^xS{x\0)e-'P-''. (29) 



S{x\0) is the inverse of the lattice Dirac operator; i.e. it is the "quark propagator" 
computed on a single gauge-field configuration and is therefore not translationally 
invariant (cf. section 4 of ref . pO] ) . It satisfies the relation 



5(x|0)=755^(0|x)75. (30) 

The amputated Green's function is obtained from eq. ( ^7| ) 

^r{l)r(2)(p)p(Tp'(T' =5' {p)paS (P)p'7 G'p(i)p(2)(p)q/3^5 S {p)f3a-S {p)sct' (31) 

where S{p) = {S{p\0)) is the Fourier transform of the translationally invariant quark 
propagator, i.e. the Fourier transform of S'(xlO), averaged over the gauge-field config- 
urations. 
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The amputated Green's function of eq. (^) is a high-rank tensor, from which a 
more manageable function of the external momenta p can be obtained by projecting 
over all the possible Dirac structures. Let us introduce a generic Dirac projector 

P^(,)f(,)^(f«"®f(2)"), (32) 

Its application on the amputated Green's functions is defined as ^: 

TvF^,,,^,,,A^,,,^,M = {t^^^^- 0t^^};^^^^ (33) 

The trace is taken over spin and colour. The projectors Pp(i)p(2) can be worked out 
analytically for the tree-level amputated Green's function 

A^SV(.) ip)f^ff = I k^"^^'"' (4^ ^ rgir ) ^ c^«^'5«'^(r2r rgi™)l (34) 



The superscript (0) denotes tree-level. From eqs. (^3|) and (E3) one finds 



±(0) ^„^ ^ ^f,(l)„^f,(2)n^^±(0) f^-^RRR'R' 

pap' a' 



TrPf„)f(,)A^(\7^(,,(p) = {T^^'>-0Ti;,';)A^i:^^,,,{p).-^.,., 



= -[Ar2(Xr f *^-^)"r^"^^'")(Tr f (2)np(2)m') 

TA^c(Tr f (i)"r(^)™f (2)"r(2)™)], (35) 

These formulae can now be used in the specific cases of interest, namely the renor- 
malization of the various parity violating and parity conserving operators (i.e. Qj^ 
of eq. ([T^) and Q^ of eq. (|l6|) with k = 1, ... ,5). We will see shortly that, for the 
renormalization conditions we will impose on the quark correlation functions of the 
four-fermion operators, it is convenient to use projectors that obey the following or- 
thogonality conditions: 

TV Pf aJ(°) = ,5,fc {i,k = l,...,5) (36) 

TV Ff :A?°^ = 5^fc (i,fc = l,...,5) 

where A^ and ^^ are the tree-level amputated Green functions of operators 
Q^ and Q^ respectively {k = 1,...,5), and P^, P^ their corresponding projectors 
{k = 1,...,5). In the following we will denote as A^ and J\^ the corresponding 
Green functions in the interacting case. With the aid of eq. ( p5| ) we find for the parity 
conserving case 



1 



P^ " +64(ivF^(P^^ - P^^) ^ 32iV.(iV2-l) (P 



55 - rpp) 



Other possible choices are admissible; see B. 

15 



,± . X ,„ „ . 1 



^5 =^32iV,(Ar2_i)(P^^ + P^^)+ 96iV,(iV2-l)P^^ 



whereas for the parity violating one 



1 

64Arc('iVe ± 1) 



?? = -7T7777^r-7T(PvA + PAy; 



^i - ^^^]v^(^(p- - P-) - i^(]v^(p- - P-) (3«) 

^ _ (2iV, ± 1) w^^^_p 

4.2 Renormalization conditions 

We now define the renormahzation procedure of the operators of interest. The general 
principle is to impose "suitable" renormalization conditions, which are satisfied by the 

renormalized (projected amputated) Green functions A^ and j7\ at a fixed scale ^ in 
the deep Euclidean region. The renormalization condition is arbitrary. A simple choice 

is to impose that the fully interacting A^ (and J\ ), at a given scale //, are equal to 
their tree level values written in eq. (|35|); see also eqs. (^). 

We will use matrix notation for simplicity: the amputated-projected Green func- 
tions A^ and ^ denote 1x5 row vectors, whereas the projectors P and P denote 
5x1 column vectors. In this notation, eqs. ( [3^ ) become 

P±.A±(o)=/ 
f± . ^±(0) ^ J (39) 

with / the 5x5 unit matrix. For the parity- violating case, the renormalized amputated 
Green's function is given by the row vector 

^^(p/m,5') = ^^'(a/x,gg):A^(ap,5o)2^(«/«,9of (40) 

where Z^ is the quark field renormalization constant and Z is the 5x5 renormalization 
matrix (the superscript T stands for transpose). Recall that Z^ is the block diagonal 
matrix of eq. (19). In the above expression we denote by Qq = g^ia) the bare coupling 
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and by g^ = g'^ip) the renormalized one. We express the bare Green function in terms 
of a "dynamics" matrix D from which the tree-level amputated Green's function is 
factored out: 

Note that since the matrix D determines the "dynamics" of the bare operators, it can 
only mix tree-level operators with the same discrete symmetries. Thus, it is also a 
block diagonal matrix with the same block structure as Z . Prom eqs. ( p9| ) and ( p] ) 
the elements of the matrix D^ are expressed in terms of the amplitudes ^ : 



V' 



^' 



(42) 



we 



We compute J\ non-perturbatively, at fixed coupling g^ and in a given gauge, over 
a configuration ensemble. We opt for the Landau gauge. From ^ and eq. 
obtain D . 

The renormalization conditions 



A (p//^,5^ 



determine the mixing matrix Z^. This is easily seen by combining eqs. ^ 
and (Bsl) to obtain 

/ 



(43) 



z^-'^v^z^^ 

from which we obtain Z in terms of the known quantities Z^/^ and D : 

-1 



'i> 



r± 



r7±2 



V 



±T 



12|) 

(44) 

(45) 



Note that since both D^ and Z^ have the same block diagonal structure, this involves 
inverting at most 2x2 matrices. 

The quark field renormalization Z^ is also determined non-perturbatively from the 
numerical simulation. A definition which respects WI's (cf. sec. 4 of ref. ||2^ for details) 

is 

1 „ dS{p)-^ 



^■0 



Instead, we have implemented 



Z^{fia) 



12 ^ %. 



1 E^=i,4 7m sin (p^a)S' (pa 



12 



4Eu=i,4 7Msin {p^,a) 



(46) 



(47) 



in order to avoid derivatives w.r.t. discrete variables. 

The computation of the parity-conserving matrix Z^ follows similar lines. Recall 
that this matrix can be cast in the form Z^[I + A ], with Z^ the block diagonal 
matrix of eq. (p^) and A^ the sparse matrix of eq. (p3|). We find it convenient to 
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obtain Z^ and A separately. As we will show in the next section, Z^ depends on 
the renormalization scale a^i^ whereas A does not. More importantly, as we will also 
show in sect. ^ the independence of A^ from o/i is due to the fact that the lattice 
subtractions can also be determined from WI's, in the spirit of ref. [^]. Thus, our 
results on A^, obtained in the RI scheme with the NP method, can in principle be 
compared to those obtained with WI's. 

For the parity-conserving case, the fully interacting renormalized amputated Green's 
function is given by the row 1x5 vector 

A^(p/^,<7') = Z;^\aii,gl)k^{av.Ql){l^^^{glf\Z^{aix,glf (48) 

We express the bare Green's function in terms of a "dynamics" matrix D^: 

A± = A±(o)l)± (49) 

From eqs. (pg|) and ( |4^ ) the elements of the matrix D are expressed in terms of the 
amplitudes A^: 

D^ = V^ ■ A± (50) 

Therefore D^ can be computed non-perturbatively by numerical simulation, at fixed 
bare coupling g^^ over a configuration ensemble and in the Landau gauge. Once it is 
computed, we determine the mixing matrix A^ and the renormalization matrix Z^ 

A. 

from the renormalization conditions 



P^^ip/lJi^g"^ 



I (51) 



To see this explicitly, we combine eqs.(|39D, (48), (^9|) and (pOj), in order to express 
eq.(^) as follows: 

VZ)±[/ + A±"]Z±^ = I (52) 

We then proceed in two steps. Let us first rewrite the above expression as 

Dti + J:D%^i = Zl\z^-\, (53) 



i=i 



Ik 



(in matrix-component notation). We then consider the special case, with indices k,l 
chosen so that for a given fixed value k, the index / is allowed to run over the cor- 
responding "lattice subtracted" values. For example, if we fix /c = 1, / runs over the 
range / = 2, 3, 4, 5; if we fix A; = 2 or A; = 3, / runs over the range Z = 1, 4, 5 etc. With 
these choices of k,l, and given the structure of matrix A (c.f. eq. (|23|) ) , the summed 
index j also runs over the same interval as I (all other contributions involve zero ma- 
trix elements of A ). Moreover, the r.h.s. vanishes for these combinations of k,l, due 
to the block-diagonal structure of the matrix Z^ (c.f. eq.(^)). Defining the column 
vector C; = D^^ for fixed k and / running as detailed above, we obtain the equation 

Ajcf = -cf (54) 



This is a linear inhomogenous system, which can be solved for A,-. Thus, the mixing 
matrix is determined, and the subtracted correlation function can be constructed. We 
can now proceed to the second step, which is the determination of the renormalization 
matrix Z^. It can be obtained exactly like in the parity violating case; cf. eqs. ( ^4[ ) 
and (H). 



4.3 Range of validity of the RI renormalization scheme 

Having completed our discussion on the determination of the renormalization constants 



and mixing coefficients of four-fermion operators, using the NP method of ref. |20|, we 
now summarize the renormalization conditions. For the parity violating operators we 
have 



Z^ ^Tr P± . J\^Z^ ^ 



For the parity conserving operators we have 



ZT^Tr P^ • A^Z^ ^ 



ZT^Tr P=^ • A*Z 



±v±T 



Z^'^Tr P^ ■ A^ 



[I- 



A±^1Z±^ 



(55) 



(56) 



with Z^ determined from eq. (p7|). In the last equation, A='= denotes the amputated 
Green's function for the lattice-subtracted parity conserving operators. These condi- 
tions constitute a renormalization scheme, called the Regularization Independent (RI) 
scheme |20, ^ (also known as the MOM scheme). The name RI is chosen so as to 
distinguish it from the MS scheme, adopted in perturbation theory, which depends on 
the detailed choice of dimensional regularization ('t Hooft-Veltman, Dimensional Re- 
duction, etc.). In the RI scheme, the renormalized operator depends on the gauge and 
on the momenta of the external states |2C , E^] . Unlike the MS, the RI scheme does not 
depend on the regularization. Since the renormalized matrix element must be multi- 
plied by the Wilson coefficient CwifJ') (cf. eq. ([l|)) the latter must also be calculated in 
the same gauge and with the same external states in order to obtain a renormalization 
group invariant result. The Wilson coefficient is known in next-to- leading order con- 
tinuum PT y]; thus, its matching with the renormalized matrix element is accurate to 
that order. Clearly, for the perturbative calculation of Cwi/-!-) to be reliable, we must 
ensure that /j, ^ Aqcd- 

There is also another reason for which this condition must be imposed: the RI 
conditions are such that the renormalized operators transform according to irreducible 
representations of the chiral algebra. This is achieved by imposing that the projections 
of the renormalized operators on the "wrong" chiral structures vanish (cf. eq.(|56D). But 
this not true when there are other causes of chiral symmetry breaking, either due to the 
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explicit presence of mass terms or due to spontaneous symmetry breaking in the chiral 
limit. Such effects, present in QCD, would modify the r.h.s. of the renormalization 
conditions ( |55|) and (|56| ) by form factors proportional to the quark mass m, the chiral 
condensate {'4)ip) and the inverse scale ;U^^. These terms die off in the large momenta 
region, ^ » Aqcd- 

On the other hand, the non-perturbative renormalization of the four-fermion oper- 
ator in the RI scheme involves computations of bare matrix elements at finite lattice 
cut-off. In order to have good control of the discretization errors, we must also ensure 



that fi <^ 0{a ^). Thus, as already pointed out in |2C], the RI scheme is applicable at 
couplings Qq and lattice sizes for which there exists a window of fi values satisfying 

Aqcd < ^ < 0{a-^) (57) 

The discussion leading to the necessity of the above window of /u values is based on 
the assumption that the operator renormalization is performed at fixed UV cutoff (i.e. 
the gauge coupling is fixed in the numerical computation). One can relax the bound 
of eq. ( |57[ ) by performing a sequence of computations at bare couplings goia), go{sa), 
go{s'^a), • ■ ■ corresponding to increasingly coarse lattice spacings (i.e. the scaling factor 
s > 1; typically s = 2). This method of lowering the renormalization scale /u in a 



controlled NP way is reminiscent of the renormalization procedure of ref. |26], which in 
turn was inspired by refs. |jl^. We outline it here in schematic fashion; wave function 
renormalization, operator subtractions etc. will be omitted for the sake of clarity. 
Realistically, with present day resources, this procedure can be iterated only two or 
three times. 

We first select a small coupling constant (^q, so that we are safely in the perturbative 
region (i.e. gQ is smaller than current values in standard QCD simulations). For such 
small lattice spacings, the window of eq. ( |57| ) is easily satisfied by both momentum 
scales /i and fi/s used in the following. We now chose a lattice of N lattice sites in each 
direction (N is only limited by current computational capabilities), in such a way that 
a/x <^ 1, in order to avoid (UV) lattice artifacts, and at the same time large enough to 
have negligible finite size effects; i.e. a///s ^ 1/A^. With small gQ ~ 0.7, and N ~ 32- 
48, this lattice is adequate for simulations of quark correlation functions deep in the 
perturbative region, but cannot accomodate hadronic quantities, due to the smallness 
of its physical volume. Once the RI scheme has been implemented at the scale n, the 
renormalized correlation function of a given operator is known at any scale p satisfying 
eq. (13): 

A (^,5'(/t^)) = Z {afi,gl (a)) A [ap,gi (a)) (58) 



In particular, we define 



AiEEA(i,5'(^)) (59) 



i.e. the renormalized vertex at a scale p = /u/s. We now increase the lattice spacing to 
a' = sa. This is done in practice by tuning the bare coupling (^q, on a bigger (coarser) 
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lattice. For simplicity, since we assumed that finite size effects are negligible both at 
^ and ^/s, we keep A'^ fixed. On the coarser lattice we compute the bare correlation 
function at momentum p' = fi/s (i.e. for a'p' = a^). The renormalization constant at 
the scale ji (for cutoff a') can then be obtained from the equation 

Ai = Z (a'/i, gl {a')) A (a'p', gl [a')) (60) 



since the l.h.s. of this expression is known from the calculation on the finer lattice 
(eq. (|59D). With the renormalization constant Z (a'/i,5Q (a')) thus obtained on the 
coarse lattice, we now compute the bare correlation function at momentum p" = /i/s^ 
(i.e. a'p" = ajj,/ s). For this momentum we then have 

A (^,ff'(M)) = Z [a'ii,gl (a')) A [a'p",gl (a')) (61) 

Thus the renormalized correlation function is now known at a lower momentum scale 

K2^kl\,g\ix)] (62) 



One can now repeat the cycle: going to a coarser lattice a" = s^a, we compute the bare 
correlation function at momentum p" = fi/ s^ (i.e. a"p" = ajj,). The renormalization 
constant at scale // is then obtained by solving 

A2 = Z (aV, gi {a")) A {a"p", gl {a")) (63) 

for Z [a" fi, gQ (a")) . The bare correlation function is then computed at momentum 
p'" = n/s^ (i.e. a"p"' = afi/s) and renormalized by Z at scale /i, giving the renormal- 
ized correlation function at scale 1/s'^: 

A (^,/(m)) = Z (aV,5o' (a")) A {a"p"',gl {a")) (64) 

At the end of the day, by keeping A^ fixed, we end up with a lattice coarse enough to 
contain the hadrons, and an operator, renormalized non-perturbatively at a low scale 

5 Mixing coefficients from WI's 



In 1 18 1 an alternative approach to the NP renormalization method has been applied. 



This is based on lattice WI's on quark states, in the spirit of [p!l| ] and [12, Q. The WI 
renormalization has been implemented in |18] for the non-perturbative evaluation of the 
mixing renormalization constants of the AS = 2 matrix element of K^—K^ oscillations. 
In our language, this is the renormalization of the parity conserving operator Qf. The 
WI's were derived in momentum space, the various chiral structures were projected 
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out (just like in [^I[ and the present work) and they were solved for the mixing lattice 
coefficients A^^ (for k = 2, . . . ,5). The multiplicative renormalization constant Z^ was 
obtained using the NP method of [^, [2^. In principle, both methods are equivalent 
(under certain conditions, which will be discussed below) and suffer from the same 
sources of systematic error. In practice, one may give more stable results than the 
other. A systematic comparison of the two methods is an interesting problem; for the 
AS* = 2 case it has been carried out in ref. []l^] . 

Here we wish to extend this method to the complete basis of operators. We will 
therefore write down the WI's which can be used for the extraction of all mixing 
coefficients (i.e. all the elements of the matrix A^ of eq. (p^)). We will also demonstrate 
explicitly that the WI method for the determination of the lattice mixing coefficients 
is equivalent to the RI method. This equivalence is true, up to discretization errors, 
provided the renormalization scale ^ of the RI scheme is large (cf. eq. (|57|)) and that 
we are working in the chiral limit or with a mass independent scheme. The reasons 



for requiring a large // for the RI scheme have been explained in subsect. 4.3. We now 
present briefly the reasons for which the implementation of the small m limit for the 
WI scheme is required. 

The WI holds for operators with the correct chiral properties, i.e. multiplicatively 
renormalizable operators transforming according to a well defined representation of the 
chiral algebra. By imposing its validity on the renormalized operators one can fix the 
mixing coefficients of the operators which stem from the chiral violation due to the 
Wilson term. This, however, is only true in the chiral limit, where chiral symmetry 
breaking is only due to the presence of the Wilson term in the lattice action. Upon 
introducing a mass term (soft symmetry breaking), the continuum WI's can only be 
recovered by the simultaneous redefinition of the T-product and the operator. This 
procedure is ambiguous, since a change in the definition of the T-product can be 
compensated by a redefinition of the operator. Thus, away from the chiral limit, the 
WI's do not uniquely define the operator (and, consequently, its mixing coefficients). 
This ambiguity is harmless if we can apply low-energy theorems of Current Algebra 
(small quark mass). This point has been discussed in ref. |jll|, and more recently in 

We will explicitly demonstrate in this section that the WI method and the RI 
renormalization scheme are equivalent methods for the determination of the lattice 
mixing coefficients, under the conditions discussed above. Moreover, we will show that 
the X-^'S' renormalization matrices Z and Z^ cannot be determined from WI's. Only 
finite ratios of x^'S renormalization constants of opposite parity are fixed by the WI's. 
Thus the RI scheme (or some other renormalization method) is necessary for their 
determination, even when the WI method is used in order to obtain the lattice mixing 
coefficients. 
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5.1 Ward Identities on quark states 

The validity of the above statements has been shown, in a general and elegant way, 
in ref. [^]. Our aim is much more specific: we will obtain useful WI's which can be 
used in practice for the determination of the mixing coefficients of all four-fermion 
operators of interest. The above general statements will consequently be proved in 
the context of these specific WI's, which can be implemented in numerical simulation. 
In what follows, for reasons of uniformity of presentation, we will continue to work 



with operators which have four distinct flavours (unlike ref. |15], where the flavour 
group used was SU{2) and the AS = 2 operator considered carried the strange and 
down physical flavours). The most economical way of obtaining useful Wis for the 
four flavour operators defined in eqs. (O) and ( [T^ ) consists in the following trick. We 
consider an "embedding" of our operators (with four distinct flavours) in a theory with 
five flavours. In other words, our flavour symmetry group is SU{5)l <8) SU{5)ji 0. A 
suitable fleld variation of the flfth flavour will then yield Wis concerning the operators 
of eqs. (^) and (|l^). Finally, since one of our aims is to propose WI's which can be 
used in simulations, we will be working with a small flnite quark mass. It is to be 
understood, throughout the rest of this section, that the chiral limit is to be taken in 
the end. 

We now establish our notation. For the purposes of this Section it superscripts will 
be dropped in what follows from operators and correlation functions (for example, from 
the Q s, G s and As). We will be using Green functions which are expectation 
values of the following operators: 

Gk{xo;Xi,X2,X3,X4)^^^S^ = iJl{xi)^i:2{x2)J3Qk{xo)Tp3{x3)^'>p4{x4)f 
Gk{xo;Xi,X2,X3,Xi)^^^S^ = ^plixi)^ij2{x2)J3Qk{xo)i^3{x3)^Tpi{Xi)f (65) 



for A; = 1, . . . , 5. The operators Qk and Qk on the r.h.s. are those of eqs. (|l^) and (|T6|); 
i.e. they are defined in terms of fermion fields with fiavours 1, • • • ,4. Their vacuum 



expectation values are the four-quark correlation functions of eq.(25), which are the 
starting point of the RI renormalization procedure. We will also be using five other 
operators (and their four-quark correlation functions) which are defined as follows: 

gk{xo;Xi,X2,X3,X4)^p^S^ = Tpl{xi)^ll)2{x2)j^ Qk{xo)4^3{x3)^ IpiiXijf 
Gk{xo;Xi,X2,X3,X4)ap^S^ = b^i'bixi)]^ '4'2{x2)p Qk{xo)'4'3(.X3)^ 'l[)4(.X4)f (66) 

where 

Qk = (v^5r(^)^2) {'tPsr^^^i^i) ± (i^^r^^^i) {i^3r^^'>i^2) (67) 

and k = 1, • • • ,5 stands for the combination of Dirac matrices corresponding to the 
definitions of eq. (16). In other words, Q^ is obtained from Q^ by substituting tpi 



^In the end of this subsection, we will comment on the relevance of our results to the physical case of 
three light quark flavours. 
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by '05. The operators Q^ are obtained from Qk, with the substitution (V'i,'0i) —>■ 
(7505i ipb)- The colour and spin indices {A, B, . . . and a, /3, . . .) will be dropped from 
now on for notational simplicity. We will also need the axial current and pseudoscalar 
density 



Pix) 


= Mxh^Mx) (68) 


We perform infinitesimal axial tra 
raising operator 

A = 


nsformations on quark fields in terms of the SU{5) 

f l\ 



(69) 


i^ / 


which induces transformations between the first and fifth flavours only: 


ipi ^ipi + iSa'y^ip^ 


V'5 ^ V'5 + iSaipi-/^ (70) 


Correspondingly, operators Q transform into their Q counterparts: 


Qi ^ Qi-i6aQi 


Q2 -^ Q2 + i6aQ2 


Q3 -^ Q3 + i6aQ3 (71) 


Qi -^ Q4 + i6aQi 


Qs 


-^ Q5 + iSaQ^ 





Whenever convenient, matrix notation will be used, as in subsect. 4.2. We define the 
1x5 row vectors 



Q = (^1,^2,^3,^4,^5) 
G = (—Gi,G2, 03,^4,^5) 

/-"S r^5 fb nb nb r'^\ 



(72) 



Note the minus sign in the element Gi of vector G; it has been introduced in its 
definition in order to take into account the minus sign of the axial variation of Qi; c.f. 
eq.(O). 

The lattice WI arising from this axial variation of the 1x5 row vector {Q) is 



5a{x) 



(^(xo;Xi,X2,X3,X4)) 



^ 5a{x 



-Q{xQ]Xi,X2,X^,Xi)) = (^(xo;Xi,X2,X3,X4) 



= 0-^ 

5S 
5a{x) 



(73) 
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where all points xo,xi,X2,X3 and X4 are kept separate. The variation of the operator 
Q receives a contribution {G) from the variation of Q and a contribution from the 
variation of the fermion field V'l • Integrating over x yields 

(G) + (^5) = - j d^x{g [V^A^(x) - 2moP{x) - X^ix)]) (74) 

In the chiral limit, the term 2mQP{x) on the r.h.s. vanishes, whereas the integrated 
total divergence of the axial current gives a non-zero surface term, due to the presence 
of Goldstone bosons. If we are not in the chiral limit, the 2m,oP(x) term is present, but 
the surface term from the current divergence vanishes upon integration. We will be 
considering the latter case, in order to mimic what is happening in the simulations (i.e. 
first we compute at small non-zero quark mass and then extrapolate to the chiral limit). 
The operator X5 arises from the variation of the chiral symmetry breaking Wilson term 



in the action. As shown in ref. [|ri| (see also |17] for a detailed discussion) it mixes, 
under renormalization with V^^^ and P. This mixing, determined by the requirement 
that on-shell matrix elements of the subtracted X^ vanish in the continuum, generates a 
finite renormalization of the axial current and a power subtraction of the quark mass. 
Thus, following ref. |]1^], in the above WI we will trade-off X^ for the renormalized 
expression [V^A^ — 2mP — X^], where X^ is the subtracted X^. What is of interest to 
us is that, besides the above renormalizations, the X5 insertion in the above correlation 
function also generates contact terms. They are found by looking at the flavour content 
and discrete symmetries of the specific correlation function. We find (up to Schwinger 
terms which will vanish under the integral): 

{X5{x)gk{xo;Xl,X2,X3,Xi)) = {Gj){5jk- Rjk)S{x- Xq) 

- {Gj)AijRik6{x-xo) (75) 

where repeated indices are summed; i,j = 1, • • • , 5. The notation for the various coef- 
ficients has been chosen with some foresight: The 5x5 matrix Rjk is a block diagonal 
matrix of the form of Zy, of eq. (p3). It will eventually turn out that R = Z^{Z~^)^; 



Jy ^-^ ^C£. \piJ/. iU WILL C VCliUUOi-LlJ L U.-L -LI WUU tlJ-O-L J l> Y 

i.e. R is the finite ratio of the parity-conserving to parity-violating, logarithmically 
diverging, renormalization matrices of xRS type. The 5x5 matrix Ajj is a sparse 
matrix of the form of eq. (p3|); it will eventually turn out to be the matrix defined in 
eq. (^). Finally note that the product [A^i?]jfc is a sparse matrix of the form of A. 
In conclusion, we have separated the contact terms on the r.h.s. of eq. ( [75[ ) into a first 
term which mixes the (Gj)'s as in the xRS scheme and a second term which generates 
lattice mixing. Note that in the above expression, we have only included the contact 
terms generated by the proximity of X5{x) to the operators Q{xo). In principle there 
are also contact terms arising from the proximity of X5(x) to the quark fields of the 
correlation Q at points j;i,---,X4. However, as shown in ref. ||4^ , these terms van- 
ish in the continuum limit. The reason is that the insertion of the operator X5 with 
fundamental fields (once both this operator and the quark fields are renormalized) is 
proportional to the lattice spacing. 
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We now combine eqs. ( [74| ) and (^), Fourier-transform the WI (with all external 
momenta set equal to p) and amputate the resulting correlation functions. We de- 
note by A(p) and ^ [p) the momentum-space amputated Green functions of (G) and 
(G^)- The Fourier-transformed Q is denoted by G{p). The resulting momentum-space 
amputated WI is 

r n r ^ 

A{p) [l + A^\r = -3\\p) + 2m / (fx{G{p)P{x)) X{{S-^{p)) (76) 

•' 1 

We require that the above WI, up to quark field renormalization, be the one valid 
in the continuum limit for renormalized correlation functions (operators) and that it 
be identical to the corresponding nominal WI (which is just the tree-level version of 
the above equation). This implies that A is indeed the matrix of lattice subtractions. 
Moreover, R can indeed be identified with the ratio of the renormalization matrices 
Z^{Z~^) , since this combination would renormalize both sides of the above equation 
(NB: the renormalization of Q{p) is identical to that of Q{p)-, since the two operators 
only differ by a relabelling of ^^5 as V'l)- 

We now show how the WI fixes both the lattice mixing matrix A and the matrix 
ratio R. We consider the first column of the above WI; its l.h.s. concerns the operator 
Qi and its lattice mixings. By projecting this WI with five projectors of eqs. (p^) and 



(38) we construct a linear inhomogenous system of five equations w.r.t. the quantities 
Rii and [A^iij/d = i^nAi^, with k = 2, •••,5. By solving this system, the lattice 
mixing of Qi and the ratio of the multiplicative renormalization constants Rn = 
Z\\j Z\\ of Q\ and Qi (which is just Q,\ with a flavour relabelling) can be determined. 
Similarly, a system of ten equations (using 10 projectors on the second and third 
columns of the WI) determines the lattice mixing of Q2 and Q3 (i.e. the A^^'s with 
J = 2, 3 and /c = 1, 4, 5) and the elements Rik with i,k = 2, 3. The case of Q4 and Q5 
is identical to that of Q2 and Qs (but involves the fourth and fifth rows of the WI). 

So far we have shown how Wis can determine the lattice mixing coefficients and 
the finite ratio of the x^S renormalization matrices. We will now show that these 
quantities are compatible to the ones obtained from the RI renormalization conditions, 
in the large momentum limit. We project both sides of eq. ([Tq) with the 5x1 column 
vector P, defined as 

P = (-Pl,P2,P3,P4,P5)^ (77) 

(the minus sign in Pi corresponds to that of —Gi in eq. (|7^). For large p, the last 
term on the r.h.s. vanishes. This is because the explicit tjiq factor implies that the 
integrated term has one less dimension that the others so that, at large momenta it 
vanishes faster by one power of p (see ref. [^31 )• Moreover, in this limit the inverse 
quark propagator behaves as S~^{p) = iSi7^p^ (with Si a scalar form factor and up 
to 0(0) terms). This means that in this limit it anticommutes with 75. Combining 
this result with the definition of Pfc we easily deduce that 

TrPfc^^^ = -TrPfcJ\^. (78) 
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(with k,j = l,---, 5). Thus, the WI becomes 



TrP-A 



/ + A^ 



Z^ = Trf- J\Z^ (79) 



This WI, given the RI renormahzation condition for the parity violating operator 

TV F • 3\Z'^ = I (80) 



implies the RI renormahzation condition of eq. (jSq ) for the parity conserving operators 
Q. 

We now comment briefly on the relation of these results to the more realistic case 
of three light flavours. The physical operators 0~ always mix with others of lower 
dimension. This introduces further complications which are beyond the scope of the 
present work. For the physical operators O"*", there are penguin- type contributions 
which, however, cancel. Thus, the resulting Wis are identical to the ones obtained 
here with five flavours. 

5.2 Scale dependence of the renormalization constants 

We will now specify which of these renormalization constants are divergent quantities in 
the UV limit and which are finite. First we give a general discussion of the functional 
dependence of the Z's (and Z^s) on the coupling, mass, renormalization point and 
cutoff. We are interested in the physics of light quark masses, so we assume fj, ^ m. 
The Z's are dimensionless quantities; thus, in principle, they could have the functional 
dependence Z((7q, a^u, am, m/iJ.), where m stands for the degenerate quark mass and /i 
is the renormalization scale (i.e. it stands for a generic choice of the four external quark 
momenta). It must be chosen so as to regularize all IR divergences, including those 
arising in the chiral limit. For the four-fermion operators we are considering in this 
work, it suffices to simply take equal momenta; i.e. pf = fJ-'^ (i = 1, . . . ,4 for the four 
external legs; fi is space-like.). Any regular dependence on am and a/i should drop out 
in the continuum limit a ^ 0, and is therefore treated in simulations as a systematic 
error due to the finiteness of the cutoff. These errors are, say, 0{am,), 0{afj,) for the 
Wilson action and 0{amgQ), Olafig^) for the tree level improved Clover action, used 
in the present paper. Moreover, the Z's cannot be singular in am. (e.g. have ln{am,) 
terms in the small mass limit m, <^ fx) because they would diverge not only in the 
continuum limit, but also in the chiral limit (m -^ 0). Similarly, regular terms in rn-//i 
are neglected (light masses), whereas singular ones must be absent (existence of chiral 
limit). In conclusion, for light quark masses, the functional dependence of the Z's is 
in general of the form Z{gQ,a^). 

Power counting suggest that their dependence on a/j, is singular (logarithmic). How- 
ever, this is only true of the x^S renormalization constants Z^ and Z. The lattice 
mixing coefficients A and the "ratio" Z~^Z^ can be determined from a system of equa- 
tions (the projected WI's discussed above). This implies that they can be expressed in 
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terms of bare lattice correlation functions, which do not depend on the renormalization 
scale. Thus their functional dependence is of the form Z{gQ); i.e. they are finite. 

Finally, we can lift the mass degeneracy, introducing small mass differences 6m. 
This could in principle introduce dependences like powers of a6m and ln{mi/mj) in 
the Z's. However, as pointed out in [^], none of these survive: The regular terms (e.g. 
powers of a6m) vanish in the continuum limit a ^ 0; i.e. they are part of the usual 
0{am) and 0{gQam) discretization errors. Singular (power or logarithmic) dependence 
is not allowed by the requirement that the chiral limit of the theory (mi — > 0) exist and 
it be well defined. Thus, lifting the mass degeneracy does not spoil the renormalization 
pattern discussed above. 

5.3 Identities between renormalization constants 

In this section we derive useful identities which relate some parity conserving renor- 
malization constants of 5 = 1 operators to renormalization constants of their S = —1 
counterparts. These identities are formally exact, but are only approximately satisfied 
in practical computations. Thus, they are useful tests of the reliability of our results. 
In order to derive them, we first need to show that, once the operator basis has been 
renormalized using the RI scheme, any other basis of renormalized operators (formed 
by linear combinations of the original renormalized operators) also satisfies the RI 
conditions. This is very straightforward. In the first basis, the RI scheme consists in 
the following properties of the projected amputated Green functions A: 

TrPA(°) = I 

A = ZA (81) 

Tr PA = Tr P(ZA) = / 

The above are valid at the renormalization scale p^ = /j,^. This shorthand notation 
should be clear to the reader. Now we can define a new basis of operators, obtained 
by a rotation R of the original basis. This implies rotated projected amputated Green 
functions A' and rotated projectors P', which must satisfy: 

A' = RA 
TrP'AW = / (82) 

Trivially, the last equation is satisfied, provided that 

P' = PR^ (83) 



It is straightforward to combine eqs.(|8l[)-(p3) to show that also A' satisfies the RI 
renormalization condition Tr P'A' = /. The implication of this property is that if 
we use the RI scheme in order to renormalize separately the operators Ofyy+AA] s-nd 
Oy^y.j^M, for example, the renormalized operator Q[vv+AA] + Qn/v+AA] satisfies the 



RI scheme; in other words it is the renormahzed operator Qf. In this section, we will 
make use of this property. We also note that in this section we will exclusively work 
with the lattice subtracted parity conserving operators Qk and the parity violating 



operators Qk- As discussed in subsecs. ^^ and 5A, these operators have good chiral 
properties (in other words, they can be thought of as operators in the X-^'S')- 

We consider the operator 0[yyi_^y], which only mixes with the operator OFgp_pg,. 
Thus, we renormalize it in the RI scheme to obtain 

0[VA-AV] = Z220[VA-AV] " ^230[sp_ps] (84) 

The same renormalization pattern is obeyed by operator Ofy^_^y^ (it only mixes with 
Otgp^PS]) since it only involves a relabelling of flavours: 

0[VA-AV] = ^220[VA-Ay] ~ ^230[sp^ps] (85) 

The last two equations can be combined into 



Q2 — 0[VA-AV] =t 0[v/A-Ay] — ^22 02 ^ ^23 23 



which implies that 

^22 = ^22 = ^22 

^2^3. = -^2~3 = ^23 (87) 

Similar expressions can be derived for the renormalization constants of operators Q3 
and the parity conserving ones in the x^S scheme. In practical simulations, we expect 
that these identities are well satisfied in the parity violating case. For the parity 
conserving one (which also involves lattice subtractions) the agreement should only be 
approximate. 

We note the above proof rests onto two crucial properties of the operators 0[y^_^y] 
and OFgp_pgi concerned: (i) they mix with each other only; (ii) they transform into 
each other under Fierz transformations in Dirac space (cf. eq.(^)). These two proper- 
ties determine the mixing pattern of eqs. (p^) and (^g). This proof cannot be extended 
to operators Q4 and Q5 . The second requirement is satisfied by the pair of operators 
0[sp+PS] arid OFgp^pg^ — O^^- see eq. (^). However, the first requirement is not 

satisfied, since operator 0[sp+p51 mixes not only with O/^p , p^, — 0|^^, , but also with 
0[sp+ps] ~ Orj^fi (cf. Table |l|). Thus, no property analogous to eq. ( ^7|) can be found 
for the operators Q4 and Q5 (or a linear combination of them). 

6 Numerical results 

In order to test the feasibility of these ideas, we have computed the renormalization 
constants, in the RI scheme, of the complete basis of operators Qk and Qk, for k = 
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l,---,5. This was done in the quenched approximation, for the Wilson and tree- 
level improved Clover action. We have performed simulations at three values of the 
gauge coupling, namely /? = 6/(7q = 6.0, 6.2, 6.4. Our results were obtained at finite 
quark masses and extrapolated to the chiral limit Kc- The specific values of all lattice 
parameters used in these simulations can be found in tab. |3[ We also give in the table 



/? 


6.0 


6.0 


6.2 


6.2 


6.4 


6.4 


Action 
# Confs 
Volume 


C 

100 
163 X 32 


W 

100 

163 X 32 


C 

180 
16^ X 32 


W 

100 

16^ X 32 


C 

60 

243 X 32 


W 

60 

243 X 32 


K 


0.1425 
0.1432 
0.1440 


0.1530 
0.1540 
0.1550 


0.14144 
0.14184 
0.14224 
0.14264 


0.1510 
0.1515 
0.1520 
0.1526 


0.1400 
0.1403 
0.1406 
0.1409 


0.1488 
0.1492 
0.1496 
0.1500 


KlQ 


0.14551 


0.15683 


0.14319 


0.15337 


0.14143 


0.15058 


a-' 


2.16(4) 


2.26(5) 


2.70(10) 


3.00(9) 


4.00(20) 


4.10(20) 



Table 3: Parameters of the runs used for the NP calculation of the renormalization constants. 
We also give the critical value of the hopping parameter Kc and the inverse lattice spacing 
a~^, as quoted in ref. [|3]. 



the value of the inverse lattice spacing a~^, obtained in ref. [43| on the same dataset. 
All statistical errors have been estimated with the jacknife method, decimating 10 
configurations at a time. 

As has been explained in sec. ^ the NP method is based on the computation of 
quark Green's functions; thus gauge-fixing has to be implemented. We have worked in 
the lattice Landau gauge, defined by minimizing the functional 



Tr 



J2{u^i^) + ul{^ 



Possible effects from Gribov copies have been ignored. However, in analogy to the study 
of the effect of Gribov ambiguities on the renormalization of two-quark operators of 



ref. [15 1, we expect them to be small. 

We now present detailed results for one representative case, namely the renormal- 
ization constants of the Clover action at /3 = 6.2, as a function of the renormalization 
scale in lattice units. All results are extrapolated in the chiral limit. In Figs. || we 
plot the renormalization constants of operators Qf and Qf. We see that both xRS 
renormalization constants Z^^ and Z^-^ are scale dependent, whereas the mixing co- 
efficients Aifc (k = 2, ... 5) become more stable with increasing ^. This is what we 
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expect from the discussion of sect. ^. In figs. ^ and ^ we show similar results for the 
renormalization constants of the operators Q^ and Q^, for A; = 2, 3, whereas in figs. |^ 
and m we show those for k = 4, 5. Most x^-S renormalization constants of the matrix 
Z^ display a marked ^-dependence due to a non-zero anomalous dimension, whereas 
the lattice subtraction coefficients of the matrix A^ are roughly scale-independent in 
the window a/i S [1,2]. Compared to the others, the renormalization constants A^, 
Ag^ and A43 have a more pronounced variation with the scale /x; nevertheless they are 
reasonably flat in the same window. 

In order to give a more quantitative flair of our results, we present the renor- 
malization matrices, for all actions and couplings, in Appendix ^. Only results at the 
renormalization scale /i ~ 2 GeV are presented. Note that the identities derived in sub- 
sec. 5^ for the x^S renormalization constants of operators Q2 ^^d Q3 (c.f. eq.(|87|)) 
are well satisfied. For the parity conserving operators Q2 and Q3 similar identities 
appear to be approximately true. At the scale /i ~ 2GeV, these results can be directly 
used in the computation of the matrix elements of the corresponding operators. A 
first implementation of these results in the calculation of various S-parameters can be 
found in refs. [^, |30| . 

The validity of our results should be confirmed independently from Wis, computed 
with the same actions and at the same couplings. Work in this direction is in progress. 
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A Notation and conventions 

The 16 Euclidean Dirac 4x4 matrices, which form a complete basis, are denoted by 

r = {I, 7a<,^a<., 7^75,75} = {5,y,r, A^}, (89) 

where 



75 = -^^<^^lupXl^Jilulpl\ = -70717273 



c^M^ = ^[Tm'Ti^] (90) 

with e^ypx the completely antisymmetric rank- four pseudotensor with £0123 = +1- The 
Euclidean Dirac matrices satisfy the following properties 

{lp,ly} = 2V' 7m = 7a., 75 = 75- (91) 
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We also define the dual "sigma" matrix 



1 



cr. 



^u 



-.^fiupXf^pX = l5(^tj.u 



T 



(92) 



The helicity projectors are, as usual, 



L, 



7/^(1-75), 7p 



Rb 



7^(1+75)- 



and 



L^{1- 75) 



^=(1+75). 



Repeated T matrices imply summation of their Lorentz indices (if any); for example 
VV = E^7m "^ 7/i, y^ = J2^,'y^l'^ 7/^75, etc. Note however, that TT = a^u ^ cr^u 
and TT = a^^ ® a^^ means summation over the 6 independent a^^ matrices (e.g. 
summation over ^ and u with fi > u). 

The colour group is SU{Nc) with Nc = 3. The Gell- Mann group generators 
are denoted by t"',a = 1,...,N^ — 1. Fermion fields ip carry spinor and colour 
indices. Latin uppercase letters denote colour indices in the fundamental representation 
(A, B, . . . = 1, . . . , Nc), whereas Greek lowercase letters stand for Dirac spinor indices 
(a, (3, . . . , p,a, . . . = 1, . . . , 4). The letters m and n are reserved for Lorentz indices 
running on the set of Dirac matrices (as in r{i)™p{2)m-j according to the following 
convention: if we are dealing with pairs of Dirac matrices such as SS, SP etc., the 
index m is absent. If we have pairs like VV, VA etc, ra runs over the four values of 
7^. In the case of the pairs TT and TT, m runs over the six independent values of the 
Dirac matrices cr^jy and ir^jy. 



B Fierz transformations in Dirac space 

In this Appendix we gather several useful formulae concerning Fierz transformations 
in Dirac space; colour indices are ignored. We express the Fierz transformation of the 
Dirac indices of a four fermion operator as follows 

r (g) n^^ 



T(^T 



(93) 



ra/3 ® T.,5 ^ [r (g) r]"^ = r^^ ® r^^ 

The Euclidean Fierz transformed Dirac tensor products [r (g) T] ^ can be reexpressed 
as a linear combination of the complete set of the original tensor products T (g) F, by 
exploiting the completeness of the set of Dirac matrices. One has 
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5. 
V ' 
T 
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P 



<v 

T 

.A 
P 



(94) 



32 



The overall minus sign is due to the anticommutativity of the Fermi fields. Our Dirac 
r matrices are normalized as 



SS = 1, VV = 4, TT 



-6, AA 



-4, PP = 1 



(95) 



with the summation over Dirac indices understood. Recall that the sum in TT runs 
only over the six independent a^^, matricesP^ 

Prom the above equation, we can easily derive the following useful identities, con- 
cerning parity conserving operators: 



[F (g) F + A (g) Af"" 
[F (g) F - A (g) A]^^ 

[S<^S-P® Pf^ 


= F(g)F + A(g)A 
= -2[S(^S-P®S] 

i 

= --[5(gS + P(gP-T 
= -2[S^S + P^P] 


[S(S)S + P(S) Pf'' 

[SO5 + POP-T0 Tf° 



(96) 



'T] 



Moreover, by substituting ijj^ — > 75'04 in the above equations, we derive, for the parity 
violating operators Qf, Q2 and Qf of our basis: 



[F g) A + A «) Vf"" 
[F g) A - A (g Vf° 


[S®P- 


- P «) Sf"" 


[5 P + P (g Sf^ 


5 (g P + P (g 5 - T (g f 



Fd 



F(g) A + A(gF 
-2 [5 (g P - P ® 5] 

-\[V®A-A®V] 



S®P+P®S 
S + P®P] 



(97) 



2 

-2[S 



T®T 



Finally, from eq. ( |94D we can easily derive the five linear combinations of operators 
which are eigenstates of the transformation with eigenvalues ±1; they are: 



[S®S + P®P + T 
\V(^V + A' 



2(5 (g 5 - P (g P) - (F (g F - A 
2(5 5 - P P) + (F F - A 




S®S + P(^P--T 





5A]^^ 
A)f- 
A)f- 

Fd 



+ [5®5 + P®P + T(gT] 

+ [F F + A (g A] 

+ [2(5 (g5-P(gP)-(F(gF-A< 

- [2(5 (^S-P(^P) + {V(^V -A^ 



S(^S + Pi 



P--T 

3 



Thus, in Dirac space, three linear combinations are Fierz eigenstates with eigenvalue 
+1 and two are Fierz eigenstates with eigenvalue —1. 



(98) 

.A)] 
.A)] 



^^ Somctinies T and A are defined with an imaginary i prefactor in order to have a positive normalization. 
This would imply a change of sign in the T and A rows and columns of the F matrix. This convention is 
not adopted here. 
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C Numerical results 

In this appendix we present the renormahzation matrices for both the Wilson and 
Cfover action at the three coupUngs considered. Ah results are in the chiral limit. For 
each P value, we show results at a scale afi (in lattice units) such that fi ~ 2GeV. 
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C.2 Wilson action, /3 = 6.0, /i^a^ = 0.964 
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C.3 Clover action, /? = 6.2, p,^a^ = 0.617 
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C.5 Clover action, /3 = 6.4, p,^a^ = 0.313 
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C.6 Wilson action, /3 = 6.4, ^^a^ = 0.313 
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Figure 1: Renormalization constants (at (3 = 6.2; Clover action) in the chiral limit as a 
function of the renormalization scale: (a) Zii (o) and Zii (□); (b) A^^, i = 2, . . . ,5. 
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Figure 2: Renorinalization constants (at (3 = 6.2; Clover action) in the chiral limit as a 
function of the renormalization scale: (a) Z^ (o) and Z22 (n); (b) Z^ (o) and ^23 ('-')5 (c) 
Z3+ (o) and Z3+ (D); (d) Z3+3 (o) and Z+ (D). 
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Figure 3: Renormalization constants (at (3 = 6.2; Clover action) in the chiral limit as a 

function of the renormalization scale: (a) A^^, i = 1,4,5; (b) A^j, i = 1,4,5. 
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Figure 4: Renorinalization constants (at (3 = 6.2; Clover action) in the chiral limit as a 
function of the renormalization scale: (a) Z^ (o) and Z44 (□); (b) Z^ (o) and Z^^ (□); (c) 
Z+ (o) and Z5+ (D); (d) Z+ (o) and Z5+ (D). 
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Figure 5: Renormalization constants (at (3 = 6.2; Clover action) in the chiral limit as a 

function of the renormalization scale: (a) Ajj, i = 1,2,3; (b) A^j, i = 1,2,3. 
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